Holographic power-law traps for the efficient production of Bose-Einstein condensates 
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We use a phase-only spatial light modulator to generate light distributions in which the intensity 
decays as a power law from a central maximum, with order ranging from 2 (parabolic) to 0.5. 
We suggest that a sequence of these can be used as a time-dependent optical dipole trap for all- 
optical production of Bose-Einstein condensates in two stages: efficient evaporative cooling in a trap 
with adjustable strength and depth, followed by an adiabatic transformation of the trap order to 
cross the BEC transition in a reversible way. Realistic experimental parameters are used to verify 
the capability of this approach in producing larger Bose-Einstein condensates than by evaporative 
cooling alone. 

PACS numbers: 37.10.Jk, 67.85. Hj, 42.40.Jv 



I. INTRODUCTION 

In recent years, several techniques have been developed 
with the aim of producing arbitrary, time-dependent op- 
tical potentials for ultracold atoms. A fast-scanning 
acousto-optic deflector can "paint" optical potentials 
with various geometries [T|-[5]. Alternatively, the diffrac- 
tion of a laser beam by a spatial light modulator (SLM) 
can produce both discrete 4 8 and continuous [9HI~2] 
holographic traps. Arbitrary trapping potentials are par- 
ticularly appealing in view of quantum information pro- 
cessing |13[ I14j and quantum simulation with neutral 
atoms [T5], the study of superfluid flow in engineered 
waveguides [16], and for matter-wave interferometry. 

To date, whenever light patterns produced with these 
techniques are used to trap a quantum gas, the experi- 
mental sequence requires an additional trap (optical or 
magnetic), in which laser-cooled atoms are evaporatively 
cooled to quantum degeneracy before being transferred 
into the optical potential of choice. This procedure would 
be significantly simplified if the laser-cooled atoms were 
trapped directly in a holographic trap and brought to 
degeneracy via a dynamic transformation of the poten- 
tial. The required apparatus is also very compact: all 
that is needed to achieve quantum degeneracy and any 
subsequent manipulation of the quantum gas in the holo- 
graphic trap are two crossed laser beams each modulated 
by an SLM. In this paper we present power-law holo- 
graphic traps suitable for this purpose, and we devise a 
sequence comprising a first stage of evaporative cooling, 
followed by a second stage in which the BEC transition 
is crossed reversibly by means of a change in trap order, 
similar to the one described in |17) . 

The paper is organised as follows: in Sec. [ll]we dis- 
cuss the method by which one can program an SLM to 
generate power-law intensity distributions, and show the 
actual light patterns produced by the physical device. To 



our knowledge, this is the first time that a linear (order 1) 
power-law optical trap, analogous to a quadrupolar mag- 
netic trapping potential, is suggested. We demonstrate 
the feasibility of using a sequence of such light patterns 
as a time-dependent trapping potential, by verifying that 
the light level does not fluctuate as the SLM is refreshed. 
In Sec. |III| we lay out the adiabaticity conditions, with 
respect to the atomic motion in the trap, which any se- 
quence must satisfy. Sec. IV details a fast and efficient 



evaporative cooling by trap deformation, while in Sec. [V] 
we analyse the adiabatic transformation of the trap in the 
vicinity of the BEC transition. We find that a large gain 
in phase-space density can be achieved during this last 
stage, with no significant atom loss. This contributes to 
the overall efficiency of the scheme in terms of the frac- 
tion of atoms left at degeneracy. 



II. GENERATING HOLOGRAPHIC TRAPS 

Using an iterative Fourier transform algorithm and an 
SLM, we can create arbitrary light patterns, including 
power-law intensity distributions / (r) of order a: 



I(r) = 




r < r , 
r > r . 



(1) 



Atoms in this light field will, via the AC-Stark effect, 
experience a conservative trapping potential for r < r$: 



U(r)=A[ — 
Jo 



(2) 
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where ro is the trap radius and A = 37rc 2 r/(0)/2o;ol , 5| 
the trap depth [18]. Here we assume light detuned by 5 
from the atomic transition frequency lo q . T is the natural 
linewidth of the atomic transition and c is the speed of 
light. 

A phase-only SLM consists of an array of liquid-crystal 
pixels which can be individually oriented. These can be 
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FIG. 1. Left: Far-field intensity distribution of a phase- 
modulated laser beam. The plane is divided into a signal 
region (inside the circle) containing a power-law distribution 
of order 2, and a noise region where the intensity is uncon- 
strained. The bright spot in the bottom right is undifiracted 
light. This intensity distribution is generated holographically 
by illuminating a phase retardation pattern (right) on a Spa- 
tial Light Modulator. 



used to manipulate the profile of a laser beam in the 
far-field by imparting a controlled, spatially-varying re- 
tardation on the light. The far-field pattern can also be 
realized by focusing the phase-modulated laser beam us- 
ing a lens. 

For a given target intensity distribution in the focal 
plane of the lens, we use the Mixed-Region Amplitude- 
Freedom (MRAF) algorithm 11 to calculate the optimal 
phase pattern in the SLM plane. This algorithm performs 
an optimization procedure on a guess phase pattern, by 
iteratively transforming the light field between the SLM 
plane and the focal plane of the lens (i.e. the output 
plane) by means of Fast Fourier Transforms (FFTs). Af- 
ter each FFT, the desired intensity of the light field is 
enforced on both planes. In particular, the MRAF algo- 
rithm divides the output plane into two regions: a signal 
region in which we restrict the intensity to match our 
target pattern, and a noise region in which the intensity 
is unconstrained. This separation allows for increased ac- 
curacy in the signal region. Phase freedom in the output 
plane is permitted, while the intensity profile in the input 
plane is that of the laser beam illuminating the SLM. The 
algorithm allows accurate and fast (a few hundred itera- 
tions) calculations of phase patterns to generate almost 
any arbitrary light pattern in the signal region of the 
output plane. An example of a phase pattern that gen- 
erates a power-law intensity distribution of order a — 2 
is shown in Fig. [I] 

The optimal retardation pattern calculated using the 
above method is then applied to an SLM (BNS P-256), 
which contains 256 x 256 24^m nematic liquid crystal 
pixels. The SLM is illuminated with 1060 nm light gen- 
erated by a diode laser, and the diffracted beam, focused 
by an / = 50 mm achromatic doublet lens, is detected 
by a CCD camera. The MRAF-optimised phase pattern 
does not always produce an accurate intensity pattern 
when applied to a physical device, due to imperfect device 
response and aberrations in the optical system [TTl IT2"] . 
However for the case of simple power-law intensity dis- 
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FIG. 2. Profiles of power-law intensity distributions with a = 
0.5 (top) and a = 2 (bottom). The measured curves (open 
circles) are extracted from images such as the one shown in 
Fig. [lj while the predicted curves (solid lines) are determined 
by the MRAF algorithm. 



tributions, we have found that the output is smooth and 
accurate. In particular we have generated power laws 
with a ranging between 2 and 0.5. Fig. [I] shows the 
light pattern in the output plane as detected by the CCD 
camera, while Fig. [2] shows the intensity profiles of a = 2 
(parabolic) and a — 0.5 distributions. At fixed order, we 
have also varied the radius ro between 27 /^m and 103 /im. 
The size of the signal region is chosen so that the uncon- 
trolled intensity in the noise region is well separated from 
the power-law pattern. This minimises the effect of the 
noise region on atoms loaded in the power-law trap. 

By applying a sequence of phase patterns to the SLM, 
dynamic optical traps can be generated. It has been 
noted previously |19j that dynamic light patterns gen- 
erated by ferroelectric liquid crystal SLMs can be sub- 
ject to substantial intensity flicker due to the changes in 
the state of individual pixels. We find that this problem 
also exists in our nematic liquid crystal SLM, but can 
be solved by careful implementation of the MRAF algo- 
rithm. To achieve dynamic power-law optical traps, we 
apply a sequence of phase patterns producing power laws 
withr a going from 0.5 to 2 in steps of 0.1 at 25 Hz. In 
the lower trace in Fig. [3j intensity flicker occurs when 
the order changes from 0.5 to 0.6 and also to and from 
order 1. This is because the initial guess phase used as 
the input to the MRAF algorithm for orders 0.5 and 1 
happens to be different from the one used for all the 
other orders. Since the MRAF-optimised phase pattern 
strongly depends on the initial guess, a significant change 
is incurred in moving from one of these orders to any of 
the others. In particular 99.7% of the SLM pixels change 
the phase shift they impart to the light, by 0.77T on av- 
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FIG. 3. Substantial nicker can occur as the SLM switches 
between patterns (lower trace). By minimising the level of 
phase change per SLM pixel, this flicker can be removed (up- 
per trace). 



erage. To overcome this, the retardation patterns are 
recalculated using the same guess phase for all orders. 
This results in a sequence in which 70% of pixels change 
their phase shift at each step, but only by 0.0257T on av- 
erage. As shown in the upper trace of Fig. [3j we can no 
longer measure any flicker between consecutive patterns. 

To ensure that the two-dimensional intensity patterns 
generated by the SLM are suitable for three-dimensional 
optical trapping of cold atoms, we verified that the beam 
maintains its power-law profile for a propagation of at 
least 300 /Ltm away from the focal plane of the lens. This 
distance is greater than the width of all the power-law 
distributions we realised, hence two perpendicular power- 
law beams, each tailored by an SLM, can provide three- 
dimensional confinement as discussed later in this paper. 
Alternatively, a light sheet can be used to add confine- 
ment along the axis of the SLM beam. 



III. ADIABATICITY CONDITIONS 

In the previous section, we showed that the size and 
order of a power-law trap can be dynamically varied by 
refreshing the phase pattern on the SLM. Moreover, the 
trap depth can be controlled by varying the optical power 
illuminating the SLM. This gives full flexibility in design- 
ing sequences. How fast a sequence can then be imple- 
mented is limited by the requirement that changes in the 
trapping potential must be slow compared to the ther- 
malisation time, which is determined by the elastic col- 
lision rate. Additionally, changes must be slow with re- 
spect to the motion of the atoms in the trap. In the 
case of a harmonic trap characterised by a frequency u>, 
any compression or expansion of the trap must satisfy 

du I u 



baticity condition can be formulated for each of the three 
parameters (radius, depth and order), given that they 
can be varied independently. The condition therefore is 
that the relative change of a parameter per unit time be 
much less than the characteristic oscillation frequency of 
an atom in the trap. For the traps considered in this 
work, this is of the order of ■J A/ m/r where m is the 
atomic mass [3TJ. Hence the adiabaticity conditions can 
be written as: 



Ak/k 1 
At ' r V rn 



(3) 



<< UJ 



Similarly, in our generic power-law trap such an adia- 



where k — A,r ai a. The SLM can vary a and r in 
small but discrete steps, and it takes At w 10 ms to 
change from one pattern to the next, as estimated by the 
duration of the flicker shown in Fig. [3] For the trans- 
formations considered below, Aq = 0.1, Ar < 3 //m, 
and the characteristic oscillation frequencies are several 
KHz. These parameters ensure that Eq. ([3| is satisfied 
for k = ro, a. The trap depth on the other hand can be 
varied continuously, hence AA can be taken as the to- 
tal change over a sequence, and At as the total duration 
of that sequence. We expect the evaporative cooling se- 
quence and the adiabatic transformation described below 
to last about one second, for which Eq. |3]) is satisfied. 



IV. EVAPORATIVE COOLING SEQUENCE 

In the remainder of the paper, we assume that two 
power law-shaped laser beams are crossed to form a 
three-dimensional trap. We assume for simplicity that 
this trap is spherically symmetric (which can be realized 
by modifying the aspect ratio of the light distributions 
shown in Sec. [IT]), but our results are also valid in the 
case of asymmetric traps as long as the trap order is the 
same in all directions. 

As for any all-optical scheme, our approach is suitable 
for all internal states, atomic mixtures and molecules. In 
the following we use 87 Rb as an example. Laser-cooled 
atoms from a magneto-optical trap can be captured into 
a broad (radius r = 50 fJ,m), 910 /iK deep parabolic 
(a — 2) optical trap. This is the deepest trap shown in 
Fig. [4^,. Following Eqs. (JlJ and ([5]), an optical power 
of 14 W in each beam is needed for such a trap at a 
wavelength of 1060 nm. 

From experiments in which similar conditions have 
been realized in practice |22) . we take as our starting 
point N — 3 x 10 5 atoms with a temperature T — 91 /iK, 
resulting in a phase-space density D = 10~ 4 and an elas- 
tic collision rate k ~ 2000 s _1 . This elastic collision rate 
is sufficiently high that a first stage of adiabatic compres- 
sion is not necessary. However this could be included (by 
reducing the trap radius and increasing the depth at con- 
stant power) if the starting conditions required it, hence 
evaporation can be optimized for a broad range of exper- 
imental parameters. We note moreover that the chosen 
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FIG. 4. The evolution of the trap profile for creating a Bose-Einstein condensate, (a) A power-law trap of order 2 is decreased 
in depth and width such that atoms are evaporatively cooled at constant elastic collision rate, (b) Following evaporation, the 
trap is adiabatically transformed from order 2 to order 0.5 to reach the critical temperature at the end of the sequence. The 
inset shows in greater detail the final stages of the adiabatic transformation. 



initial atom number is conservative, and that increasing 
the trap size results in more atoms being captured [22], 
ultimately leading to a larger condensate. 

A possible evaporation sequence is shown in Fig. |4^. 
The trap depth is gradually lowered to force evaporation, 
while the trap radius is adjusted to maintain optimal con- 
ditions. In particular, we choose to keep the elastic colli- 
sion rate constant, as opposed to selecting a sequence of 
runaway evaporation. Given our starting conditions, this 
ensures that the atomic number density does not increase 
to the point where three-body recombination losses be- 
come significant. The trap order is kept at 2 because 
smaller values would result in decreased efficiency of the 
evaporation process [23) . The efficiency is defined as: 



7 



d\nD 
din A 7 ' 



(4) 



which is equivalent to N oc 1/Dy oc T^V^ /Ny , where 
V is the effective volume occupied by the atoms. This 
relation can be used to express the dependence of N on V 
and T, which is then substituted in the condition for the 
clastic collision rate n oc T^ 2 N/V = constant. Next we 
express T and V in terms of the trap depth and radius re- 
spectively: T oc A and V oc 7q, assuming the truncation 
parameter r\ — A/(ksT) is kept fixed throughout evapo- 
ration. The condition for constant k finally becomes: 



Ai 



-2 



.37 



= constant, 



(5) 



which is used to plot the intermediate steps of the evap- 
oration sequence in Fig. [4^,. For this we assume 7 = 3 



and 77 = 10 (i.e. evaporation near stagnation), as typi- 
cal for evaporation in optical traps [22, 24 26 . Similarly 
to these experiments, we expect evaporation to work on 
timescales of seconds. Hence the SLM needs to be re- 
freshed at a rate less than 20 Hz to produce the inter- 
mediate steps. Given that this is much smaller than the 
characteristic trap frequency, we do not expect significant 
parametric heating 27 from residual intensity flicker at 
the refresh rate. 

At the end of the evaporation sequence, the final trap 
radius and depth are 20 /^m and 84 /iK respectively, at 
which point D = 0.011 and N = 6 x 10 4 atoms remain in 
the trap. The adiabatic transformation described in the 
next section is then used to bring the sample to degener- 
acy. 



V. ADIABATIC TRANSFORMATIONS 

As shown in [T7], the phase-space density of a colli- 
sional gas can be increased adiabatically and revcrsibly, 
without loss of atoms, by reducing the trap order. This 
is different from compressing or expanding a trap while 
keeping the same order, in which case the temperature 
and the number density change so as to keep the phase- 
space density constant. In the following, we consider an 
adiabatic transformation in which the trap order goes 
from a = 2 to 0.5, starting from the parabolic trap at 
the end of the evaporation sequence (i.e. the smallest 
trap in Fig. [4^,) . We define our transformation such that 
at a = 0.5 the atoms reach the critical temperature for 
Bose-Einstein condensation, i.e. T — T c and D — 2.612. 



FIG. 5. The entropy of a Bose gas in a power-law trap of order 
a as a function of (left) condensate fraction and (right) tem- 
perature above and below T c . The two dotted lines represent 
transformations at constant entropy. By transforming from 
a power-law trap of order 2 containing a gas at T/T c — 4.7 
to a power-law trap of order 0.5, the gas reaches the critical 
temperature (upper dotted line). Starting from T = T c in the 
a = 2 trap, an adiabatic transformation to a = 0.5 leads to a 
condensate fraction of about 0.6 (lower dotted line). 



For such a transformation to be reversible, it must be 
slow compared to the thermalization time. We estimate 
that the elastic collision rate during our sequence does 
not drop significantly from the value of 2000 s _1 during 
evaporation, so a sequence lasting about one second (as 
suggested in Sec. Ill I satisfies this requirement. We can 




therefore impose that the entropy S is constant during 
the transformation, with S defined by [17] : 



1.6 1.4 1.2 1.0 

Power-law a 

FIG. 6. The temperature and critical temperature decrease 
by changing the power-law order a, until they are equal. The 
phase-space density of the gas increases by about a factor 230 
to 2.612. 



lOfcsT(a) to avoid further evaporation during the trans- 
formation. Secondly, the end point of the sequence is de- 
termined by imposing a peak density n — 1.8x 10 14 cm -3 , 
which gives a moderate rate of atom loss due to three- 
body recombination of Is -1 [29]. We then use the 
critical condition nXi B = 2.612 (with the thermal de 
Broglie wavelength XdB = h/y / 2nmkBT) to calculate 
T c = 582 nK in the a — 0.5 trap. To find the trap size 
r (0.5), we consider the expression for T c in a generic 
power-law trap |30j : 



95 ,3 (z) 
9 ^ ™ 



Nkf 




T>T C , 



T<T Cl 



(6) 



where z = exp {^,/ksT) is the fugacity with \x the chemi- 
cal potential, and g K (z) — Y^jLi z " 'j K - Below T c , fi = 



and therefore z 
solving 



1. Above T c , we find z by numerically 



3| + 3 (z) = gs + 3 (1) 



(7) 



This results in the plots of entropy per particle versus 
T/T c shown in Fig. [5] for several trap orders. For the 
conservation of entropy, the value S/Nks = 8.47 at the 
critical temperature for a = 0.5 (see upper dotted line 
in Fig. [5]) must also be the entropy per particle at the 
beginning of the adiabatic transformation, i.e. for a = 2. 
This corresponds to T/T c = 4.7 and D 1 = 0.011 in the 
parabolic trap, as achieved at the end of the evaporation. 
From the conservation of entropy and Eq. |7 ) we can then 
determine z{a) and Sr (a) for the intermediate steps of 
the adiabatic sequence. 

To calculate the trap depth A{a) and size r {a) re- 
quired for this sequence, we start by imposing A{a) = 



knTr = 



{2m)- 



1 



c( 



(8) 



and we solve it for tq to find ro(0.5) = 59 fim. Having 
fixed the trap parameters at the start and the end of the 
sequence, we assume for simplicity a linear interpolation 
for r (a). We use Eq. Q, in conjunction with |r(a:) 
from entropy conservation, to determine T{a) and T c {a) 
individually, as shown in Fig. [6j The resulting trap pro- 
files are shown in Fig. |4Jd. 

The optical power needed for the final trap {a = 0.5) 
is 36 mW per beam, which is 0.26% of the power at the 
start of evaporation. Such a large dynamic range can 
be achieved by varying the power illuminating the SLM, 
e.g. with a motorized rotating waveplatc followed by a 
Glan- Taylor polarizer, and by changing the power emit- 
ted by the laser [22]. Towards the end of the adiabatic 
transformation, the trap is so shallow that the effect of 
gravity cannot be neglected. It is possible to compen- 
sate for it either with a magnetic field gradient, or by 
designing "tilted" intensity patterns. 

We see from Fig. [6] that the adiabatic transforma- 
tion achieves a significant gain in phase-space density. 
Intuitively, the gas undergoes an expansion and its tem- 
perature decreases. However the change in trap order 
causes an increase in peak density which is especially 
pronounced at the end of the sequence. We estimate 
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that only 5% of atoms are lost due to three-body re- 
combination for a sequence that lasts Is, which justifies 
our assumption of constant atom number throughout the 
adiabatic transformation. In conclusion, about 20% of 
the laser-cooled atoms loaded into the initial optical trap 
remain when the sequence terminates at the BEC tran- 
sition, which is higher than in any other all-optical tech- 
nique implemented so far. 

If evaporation is continued to the BEC transition at the 
same efficiency assumed in Sec. IV 10 4 atoms remain in 
the trap - which is only 3% of the initial atom number. 
In this case however the adiabatic change of trap order 
(see lower dotted line in Fig. [5]) provides a new method 
to cross the BEC transition in a reversible way, similarly 
to the dimple technique [3T]. From the conservation of 
entropy, we find that varying a from 2 to 0.5 leads to 
T/T c ~ 0.89. A condensate fraction N BEC /N ~ 0.6 is 
found using N BEC /N = I- (T/T c ) 3/2+3/a [T7], as shown 
in the left side of Fig. [5] hence Nbec — 6000 at the 
end of the adiabatic transformation. Our sequence could 
then be combined with in-situ imaging |32j to provide a 
reversible method to investigate the BEC transition. 



VI. CONCLUSIONS 

We have used an SLM to holographically generate 
power-law intensity patterns of different orders and sizes, 
and we have shown how a sequence of these can be used 
as a dynamic optical trap for fast and efficient produc- 
tion of Bose-Einstein condensates. Starting from realis- 
tic assumptions, we have calculated the trap parameters 
throughout the sequence. We have also presented the adi- 
abaticity criteria for a generic transformation that goes 
beyond the standard case of the compression and expan- 
sion in a harmonic trap. A future step will be the study 
of evaporative cooling and adiabatic transformations in 
more complex trap geometries. 
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